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ABSTRACT Two-component systems (TCS) comprise histidine kinases and their cognate response regulators and allow bacteria 
to sense and respond to a wide variety of signals. Histidine kinases (HKs) phosphorylate and dephosphorylate their cognate re- 
sponse regulators (RRs) in response to stimuli. In general, these reactions appear to be highly specific and require an appropriate 
association between the HK and RR proteins. The Myxococcus xanthus genome encodes one of the largest repertoires of signal- 
ing proteins in bacteria (685 open reading frames [ORFs]), including at least 127 HKs and at least 143 RRs. Of these, 27 are bona 
fide NtrC-family response regulators, 21 of which are encoded adjacent to their predicted cognate kinases. Using system-wide 
profiling methods, we determined that the HK-NtrC RR pairs display a kinetic preference during both phosphotransfer and 
phosphatase functions, thereby defining cognate signaling systems in M. xanthus. Isothermal titration calorimetry measure- 
ments indicated that cognate HK-RR pairs interact with dissociation constants (K^) of approximately 1 fiM, while noncognate 
pairs had no measurable binding. Lastly, a chimera generated between the histidine kinase, CrdS, and HKl 190 revealed that resi- 
dues conferring phosphotransfer and phosphatase specificity dictate binding affinity, thereby estabUshing discrete protein- 
protein interactions which prevent cross talk. The data indicate that binding affinity is a critical parameter governing system- 
wide signaling fidelity for bacterial signal transduction proteins. 

IMPORTANCE Using in vitro phosphotransfer and phosphatase profiling assays and isothermal titration calorimetry, we have 
taken a system-wide approach to demonstrate specificity for a family of two-component signaling proteins in Myxococcus xan- 
thus. Our results demonstrate that previously identified specificity residues dictate binding affinity and that phosphatase speci- 
ficity follows phosphotransfer specificity for cognate HK-RR pairs. The data indicate that preferential binding affinity is the basis 
for signaling fidelity in bacterial two-component systems. 
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In response to nutrient limitation, Myxococcus xanthus under- 
goes a complex developmental process culminating in the for- 
mation of multicellular structures termed fruiting bodies. As the 
fruiting body matures, individual cells differentiate into metabol- 
ically dormant and stress-resistant spores. The developmental 
process is characterized by significant changes in gene expression 
(1). Many developmentaUy regulated genes are ct^* dependent, 
thus requiring the alternative sigma factor RpoN (2). Transcrip- 
tion from cr^''-dependent promoters requires an upstream activa- 
tor homologous to bacterial NtrC proteins. Mutations inactivat- 
ing most of the 27 NtrC homologs in M. xanthus produce defects 
in development, leading to the hypothesis that NtrC-like activa- 
tors (NLAs) are involved in complex signaling networks (2-12). 
NtrC homologs are response regulators (RRs) for two-component 
signal transduction systems (TCS), which require activation by a 
histidine kinase (HK). While NtrC homologs are regulated in 
M. xanthus (2), there is a significant lack of data regarding the 
associated HKs. Indeed, M. xanthus NLAs regulate development 
through a complex cascade of gene expression (2) and may inter- 
act biochemically in functionally complex signaling pathways. 



Many organisms utilize complex signal transduction cascades 
to regulate critical cellular processes. Complex regulatory path- 
ways allow integration of multiple signals and checkpoints for the 
control of important cellular decisions. For instance, Bacillus sub- 
tilis employs multiple kinases within a complex signaling cascade 
to regulate phosphorylation of SpoOAto control sporulation (13). 
Alternatively, Escherichia coli Che A regulates the activity of both 
CheY and the negative feedback controller, CheB, to control che- 
motaxis (14, 15). In M. xanthus, the CrdS-CrdA TCS interacts 
with the Che3 chemosensory system to regulate timing of aggre- 
gation and fruiting body development (5, 6, 16). These examples 
highlight the variability in signal transduction system architecture 
and suggest that further analysis of M. xanthus two-component 
systems may provide additional insights into the regulation of 
complex signaling pathways. 

Upon sensing an appropriate ligand, a TCS HK autophospho- 
rylates on a conserved histidine. The phosphoryl group is then 
transferred to a conserved aspartate within the RR, thus regulating 
activity of the system. Many HKs are bifunctional enzymes pos- 
sessing the ability to both phosphorylate and subsequently de- 
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phosphorylate their cognate RRs. Dephosphorylation of a phos- 
phorylated RR provides a mechanism to eliminate cross talk from 
other sources which would otherwise lead to fitness disadvantages. 
Specific residues within the HK are known to be required for 
proper phosphatase activity (16-32), and specificity residues are 
known to prevent cross talk (27, 33-36). One aspect of signaling 
fidelity that has yet to be determined conclusively is whether spec- 
ificity residues dictate binding affinity to ensure cognate HK-RR 
interactions. Some examples have implicated binding or compe- 
tition as a factor in the prevention of cross talk (27, 30, 36-43) , but 
information regarding a direct link between specificity residues 
and binding affinity is limited. 

Studies with E. coli and Caulobacter crescentus established the 
basis for signaling fidelity regarding TCS phosphotransfer (44, 
45). However, by comparison to M. xanthus, those organisms 
have fewer and relatively distant TCS homologs. M. xanthus en- 
codes -127 TCS, including those utilizing NtrC homologs which 
likely arose through duplication (46), suggesting that signaling 
fidelity and specificity may be more complex in this case (43). To 
investigate fidelity and binding affinity, we performed phospho- 
transfer profiling with each of the 27 predicted NtrC RR homologs 
in combination with a subset of corresponding HKs. We also pro- 
filed phosphatase activity for ten putative HK-RR systems and 
extended our analysis to binding affinity for a subset of these pro- 
teins. Lastly we generated a chimera that switched its binding af- 
finity as well as its preferred target for phosphotransfer and phos- 
phatase specificity to its predicted new RR partner. Our results 
demonstrate a clear kinetic preference for phosphotransfer and 
phosphatase activities and correspond directly with binding affin- 
ity. The results support the model that binding affinity is con- 
trolled by specificity residues which thereby dictate kinetic prefer- 
ence. 

RESULTS 

Bioinformatic analysis of signaling systems utilizing NtrC ho- 
mologs in M. xanthus. The M. xanthus genome encodes a mini- 
mum of 127 HKs and 143 RRs, 27 of which have domain topology 
identical to that of the E. coli protein NtrC. These proteins were 
identified using the Microbial Signal Transduction Database by 
searching for proteins with NtrC-like domain topology: each ho- 
molog possesses an amino-terminal receiver domain (REC), a 
central cr"* activation domain, and a carboxyl-terminal helix- 
turn-helix DNA-binding domain (47). In addition, all of the as- 
sociated REC domains have conserved residues required for phos- 
phorylation by a histidine kinase, indicating that these proteins 
are likely involved in signaling (48). The 27 NtrC homologs are 
distributed throughout the M. xanthus chromosome and display 
high sequence similarity, which is indicative of recent duplication 
events rather than acquisition via horizontal transfer (46). A sub- 
set of NtrC-like activators have also been annotated and studied 
previously (12). Many NtrC homologs have been studied in 
M. xanthus, and most are reported to affect development (2-12) 
(Table 1). Based on these studies, it has been suggested that NtrC 
regulators maybe involved in cross talk or cross-regulated systems 
(2, 6). 

To determine if signaling fidelity is maintained for this set of 
TCS RR proteins, as described previously (33, 36), or if there 
might be some level of cross-regulation for these systems, we as- 
sessed the phosphorylation kinetics for several HK-RR systems in 
vitro. First, we sought to identify cognate kinases for each NtrC 



homolog, as those proteins have not been extensively character- 
ized. The gene neighborhood for each RR was analyzed for in- 
stances where an HK was encoded within a putative operon. Of 
the 27 M. xanthus NtrC homologs that met our criteria, we were 
able to identify 19 HKs encoded adjacent to or within the same 
putative operon or gene cluster (Table 1), 6 of which have been 
studied previously: FrgBC (49), HsfBA (10), MrpAB (50), PilSR 
(11), Nla6S (HK)-Nla6 (RR) (12, 51), and Nla28S (HK)-Nla28 
(RR) (12, 52). The NtrC homologs CrdA and SasR are orphans but 
have experimentally determined cognate HKs, CrdS and SasS, re- 
spectively (6, 7). Six additional NtrC homologs are orphans 
(NtrC0172, NtrC3555, NtrC4240, NtrC4977, ActB, and Nla4 [3, 
8, 12]) and might be secondary targets for prototypical TCS or 
targets of orphan HKs. Genetic analyses were performed for 7 
additional RRs (NtrC4196, NtrC4261, NtrC7143, Nla7, Nlal9, 
Nla23, and Nla24 [3, 12]) but did not assess a role for the associ- 
ated putative HK. The 6 remaining RRs and their putative associ- 
ated HKs have not been assayed for function in vivo. Moreover, 
none of the aforementioned kinases, regulators, or predicted cog- 
nate pairs have been tested for specificity via phosphotransfer pro- 
filing. The predicted or experimentally determined cognate HK 
and RR homolog and known phenotypes for corresponding mu- 
tants are listed in Table 1. 

NtrC homologs are preferentially phosphorylated by cog- 
nate histidine kinases. Previous studies demonstrated that spec- 
ificity during in vitro phosphotransfer profiling reflects in vivo 
phosphotransfer preference for a given HK-RR signaling system 
(44, 53 ). To characterize the phosphotransfer preference for a sub- 
set of M. xanthus HK-RR pairs, we first generated constructs to 
express the signaling domains of the 27 NtrC RRs and the 21 
associated HKs. Each construct was designed to express a soluble 
REC domain for each RR or a truncated form of the kinase (in- 
cluding the dimerization and histidine phosphotransfer [DHp] 
and the catalytic and ATP-binding [CA] domains [54] ) where the 
input domain was replaced by a 6-His affinity tag. Using this 
method, we were able to purify 25 of 27 NtrC REC domains and 
soluble, active forms of 10 kinases (Fig. 1): HK0938, HK1190, 
HK3812, HK4786, HK5778, HK7439, FrgB, HsfB, CrdS, and SasS. 
Activity for each of these proteins was first determined by in vitro 
autophosphorylation. For REC domains, the ability to incorpo- 
rate radioactive phosphoryl groups was determined after incuba- 
tion with [''^P]acetyl-phosphate (AcP) (see below). The activity 
for each HK was determined by its ability to incorporate radioac- 
tive phosphoryl groups when incubated with [y-^^P] ATP. The 10 
HK proteins displayed variation in both the amount and rate of 
autophosphorylation (Fig. 1; also, see Fig. SI in the supplemental 
material). For instance, HsfB reached maximal phosphorylation 
in approximately 5 min, while HK4786 did not attain peak levels 
even after 60 min. Additionally, the apparent maxima varied sig- 
nificantly between proteins: HK3812 was able to generate only 
0.38 pmol of phosphorylated protein, while FrgB generated 
2.4 pmol at 60 min. A similar level of variability has been docu- 
mented for HKs purified from E. coli, where both the rate and level 
of phosphorylation varied significantly (55), and this may reflect 
variation in the fraction of active or properly folded proteins. 

To test kinetic preference for each HK, we performed a series of 
phosphotransfer profiling experiments similar to those described 
previously (44). In these experiments, each individual HK was 
autophosphorylated prior to incubation with each individual RR. 
As a baseline, we performed phosphotransfer profiling with CrdS, 
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TABLE 1 HK-NtrC pairs in M xanthus' 



HK 

(Mxan no.) 


HK designation 


NtrC (Mxan no.) 


NtrC designation 


ntrC (nla) mutant 
phenotype 


Reference for HK 
phosphorylation or 
phosphotransfer 






0172'' 


Nti-C0172 


Development (3) 




0938 


HK0938 


0937 


NIa7 


(12) 


This study 


1077 




1078 


NIal9 


Development (12) 




1129 


FrgB 


1128 


FrgC 


Development (49) 


This study 


1166 


Nla28S 


1167 


Nla28 


Development (12) 


52 


1190 


HK1190 


1189 


NtrC 1189 




This study; 6, 16 


1249 


SasS 


1245 


SasR 


Development (7) 


This study 






2516'' 


Nla4 


Development (12) 








3214'' 


ActB 


Development (8) 




3419 




3418 


NtrC3418 










3555'' 


NtrC3555 


(3) 




3812 


HK3812 


3811 


NtrC3811 




This study 


4043 


Nla6S 


4042 


NIa6 


Development (12) 


51 


4197 




4196 


NtrC4196 


Development (3) 








4240'' 


NtrC4240 






4251 




4252'' 








4262 




4261 


NtrC4261 


Development (3) 


6, 16 


4579 




4580' 








4786 


HK4786 


4785 


NtrC4785 




This study 






4977'' 


NtrC4977 


(3) 




5123 


MrpA 


5124 


MrpB 


Development (50) 




5184 


CrdS 


5153 


CrdA 


Development (5, 6) 


This study; 6, 16 


5365 


HsfB 


5364 


HsfA 


Heat shock (10) 


This study; 10 


5778 


HK5778 


5777 


Nla23 


Development (12) 


This study 


5785 


PilS 


5784 


PilR 


Motility (11) 




7142 




7143 


NtrC7143 


Development (3) 




7439 


HK7439 


7440 


NIa24 


Development (12) 


This study 



" shown is the list of M. xanthus NtrC homologs and their predicted cognate partner histidine kinases (HK). The locus tags for each HK or NtrC RR homolog are shown. 
Alternative names for some of the homologs and phenotypes for mutants are also listed. Biochemical determination of in vitro phosphorylation of the HK is listed in the last 
column. Mxan, the locus identifier for eachM. xanthus gene in the NCBI database. 
^ Orphan ntrC gene (lacking a corresponding kinase). 
The NtrC homolog was not tested. 



for which both genetic and biochemical evidence indicates that 
CrdA is the cognate RR in vivo (6, 16). When we allowed the 
experiment to proceed for 10 s, we saw phosphotransfer from 
CrdS to CrdA (Fig. 1). However, when we allowed the same ex- 
periment to proceed for an hour, CrdS was able to phosphorylate 
three other NtrC homologs, NtrCO 1 72, NtrC 1189, and the orphan 
NtrC3555 (see Fig. S2 in the supplemental material). It is likely 
that CrdS phosphorylates these alternative targets only under the 
conditions of this assay and is not physiologically relevant. This 
result is similar to other phosphotransfer profiling experiments 
where the E. coli kinases EnvZ, CheA, and CpxA phosphorylate 
their cognates at short time points but allow transfer to alternative 
targets over longer time courses (44). Because longer time points 
are generally permissive for low-level phosphorylation of nonspe- 
cific targets, we performed all subsequent phosphotransfer profil- 
ing experiments using a 10-s time point. 

In addition to CrdS, we performed phosphotransfer profiling 
experiments using the M. xanthus kinases HK0938, HK1190, 
HK3812, HK4786, HK5778, HK7439, FrgB, HsfB, and SasS. The 
results for the phosphotransfer profiles were similar to those ob- 
tained with CrdS. Each HK phosphorylated its predicted cognate 
RR within 10 s (Fig. 1), although variability was apparent in this 
experiment. For instance, there was no observable loss of phos- 
phorylated HK5778 (HK5778~P) which corresponded with only 
modest production of the phosphorylated target, Nla23~P. In 
stark contrast, HK1190 exhibited significantly higher levels of 



phosphotransfer to NtrCl 189, as indicated by the complete turn- 
over of HKl 190~P which corresponded with production of high 
levels of phosphorylated target protein, NtrC1189~P. Collec- 
tively, the results show that each HK displays a kinetic preference 
for phosphotransfer to its predicted cognate RR and that these 
systems comprise discrete signaling pathways: HK0938-Nla7, 
HK1190-NtrC1189, HK3812-NtrC3811, HK4786-NtrC4785, 
HK5778-Nla23, HK7439-Nla24, CrdSA, FrgBC, HsfBA, and 
SasSR. 

Phosphatase profiling demonstrates specificity for bifunc- 
tional histidine kinases. Many HKs are known to be bifunctional 
enzymes capable of both phosphorylation and dephosphorylation 
of the cognate RR (17-20, 23-26, 30-32). Phosphatase activity is 
critical to signaling fidelity, as it provides a mechanism to negate 
cross talk or aberrant phosphorylation by other sources (27, 33- 
36, 43). Thus, we sought to test the fidelity of HK phosphatase 
activity by conducting a corollary assay termed phosphatase pro- 
filing. In this assay, we generated RR~P using the high energy 
phosphodonor acetyl-phosphate (AcP). We selected the 10 NtrC 
homologs that were specifically phosphorylated by their cognate 
HK in our phosphotransfer profiling assays (Fig. 1): NtrC1189, 
NtrC381 1, NtrC4785, Nla7, Nla23, Nla24, CrdA, FrgC, HsfA, and 
SasR. Each RR was labeled using [^^P] AcP as described previously 
(6, 36, 56). Each RR~P incubated alone in buffer was used to 
normalize the data for subsequent reactions (Fig. 2). Each RR~P 
was incubated with an equimolar amount of each of the ten puri- 
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FIG 1 Phosphotransfer profiling indicates specificity for cognate TCS. Each HK was allowed to autophosphorylate to maximum levels before incubation with 
an equimolar amount of each individual NtrC RR homolog. All phosphotransfer reactions were stopped after 10 s, and products were resolved using SDS-PAGE 
and analyzed on a phosphorimager as described in Materials and Methods. Each phosphorylated NtrC homolog is indicated by an asterisk. Three separate gels 
were run and processed simultaneously for each HK to enable profiling for all 25 RRs. The resulting images were assembled into rows as shown. 



fied kinases. A 5-min incubation was chosen because the assay 
requires measuring the disappearance of label, and that time point 
had been used previously to demonstrate phosphatase activity for 
several other kinases (16). 



Initially, we performed the phosphatase profiHng experiment 
using the previously described NtrC homolog, CrdA. In this ex- 
periment, CrdA was dephosphorylated only by its cognate kinase, 
CrdS. In the presence of CrdS, CrdA~P levels were reduced to 
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FIG 2 Phosphatase profihng reveals specificity within cognate TCS. Phosphatase assays were performed using acetyl-phosphate (AcP)-labeled response 
regulators as described in Materials and Methods. The amount of phosphorylated RR incubated in buffer alone was set to 100%. A decrease in the amount of 
phosphorylated RR is indicative of phosphatase activity. 



1.6% relative to the buffer control (Fig. 2). When CrdA~P was 
incubated with each of the other HKs, there was no statistically 
significant diminution of CrdA~P, indicating that CrdS phospha- 
tase activity is specific for CrdA~P. Using this phosphatase profil- 
ing assay, we observed a similar trend for each of the other nine 
NtrC homologs. Each RR~P was dephosphorylated only when 
incubated with its predicted cognate HK (statistical significance 
determined by Student's t test; P < 0.05). Four HKs (HK5778, 
CrdS, FrgB, and SasS) were able to dephosphorylate their cognate 
RR~P by 90%. The remaining 6 HKs (HK0938, HKl 1 90, HK38 12, 
HK4786, HK7439, and HsfB) dephosphorylated their targets by 
-60 to 90% (Fig. 2). The variability for phosphatase activity ob- 
served here is similar to that observed for phosphotransfer activity 
exhibited by each HK. Importantly, each HK displayed specificity 
for the same RR homolog for both phosphatase and phospho- 
transfer activities (Fig. 1 and 2). Together, these data allow us to 
conclude that phosphatase activity, like phosphotransfer activity, 
is highly specific and can be used to define cognate signaling sys- 
tems. 

A combined assay illustrates differences in phosphorylation 
kinetics for TCS. After determining the phosphotransfer and 
phosphatase specificity independently for each HK via profiling, 
we then analyzed the combined phosphotransfer and phosphatase 
functions over time in one assay. In essence, the assay measures 
the flux of phosphoryl groups through the HK-RR system. We 
tested the 10 HK-RR pairings identified above in the profiling 
experiments. CrdS~P displayed rapid phosphotransfer to CrdA. 
Subsequently, CrdS dephosphorylated CrdA~P (Fig. 3), resulting 
in a complete loss of radiolabel from the reaction by 1 min in this 
assay. Because the half-life of CrdA~P is about 1 h (6), the disap- 
pearance of both CrdS~P and CrdA~P indicates fast phospho- 
transfer and subsequent dephosphorylation by CrdS. HK0938, 



HKl 190, HK4786, HK7439, FrgB, and HsfB demonstrated kinet- 
ics similar to that of CrdS, where < 10% of the HK~P remained at 
the 1-min time point. Additionally, these HKs were able to de- 
phosphorylate their cognate RRs to various degrees over the 
5-min time course. In contrast, HK3812, HK5778, and SasS dis- 
played lower phosphotransfer rates, with more than 25% of HK~P 
being present at the 1-min time point. HK5778 displayed the 
slowest turnover, with more than 30% of HK5778~P remaining 
after 5 min. Overall, the results are consistent with the phospho- 
transfer and phosphatase profiling data yet reveal differences in 
rates that may reflect physiologically relevant aspects of TCS con- 
trol in M. xanthus. 

Cognate systems display preferential binding as measured 
by isothermal titration calorimetry. Previous work led to the hy- 
pothesis that binding affinity facilitates preferential interactions 
between TCS proteins to maintain signaling fidelity (27, 30, 36- 
43). In this study, we used isothermal titration calorimetry (ITC) 
to determine binding affinities of four pairs of TCS proteins and 
tested all cognate and noncognate pairwise interactions. We chose 
M. xanthus HK1190-NtrC1189, CrdS-CrdA, FrgB-FrgC, and 
T. maritima HK853-RR468. We chose HKl 190 because it is most 
similar to CrdS within the al helbc of the DHp domain (see be- 
low). Additionally, RR NtrC 1 189 was phosphorylated, albeit inef- 
ficiently, by CrdS at the 1-h time point in the phosphotransfer 
profiling experiment (see Fig. S2 in the supplemental material), 
suggesting the possibility of cross talk between these two systems. 
We chose FrgBC because genetic analyses confirmed its likely 
function as a cognate TCS affecting motility in M. xanthus (49). 
Lastly, we chose the well-studied yet evolutionarily distant 
HK853-RR468 TCS from T. maritima to serve as a control, and we 
determined that HK853 exhibited no capacity to phosphorylate 
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FIG 3 Time course assays for phosphotransfer activity. Cognate HK and RR pairs identified by tlie above profiling experiments were assayed in a time course 
for phospliotransfer activity. Each HK protein (black) was allowed to reach maximal phosphorylation and then mixed with equimolar amounts of cognate NtrC 
RR protein (red). Samples were collected at 0, 0.25 0.5, 1 and 5 min, resolved by electrophoresis, imaged, and quantified as described in Materials and Methods. 
The quantified results are presented as pmol of phosphorylated proteins. 



any of the 25 M. xanthus NtrC RR homologs in this study (see 
Fig. S3 in the supplemental material). 

Background-corrected, integrated thermograms for three rep- 
resentative binding reactions are shown in Fig. 4. Calorimetric 
measurements and subsequent fitting of binding isotherms al- 
lowed us to determine the affinity for each pairwise protein- 
protein interaction. For the cognate pair, CrdS-CrdA, the reaction 



curve clearly indicates an interaction with an estimated dissocia- 
tion constant (K^) of 1.4 ± 0.3 /liM. Similar results were obtained 
for each pairwise interaction of cognate partner proteins, where 
the resulting was estimated to be about 1 (Table 2). It is 
worth noting that, under the conditions of the assay, the proteins 
are unlikely to be phosphorylated. Nevertheless, experiments per- 
formed using noncognate proteins pairs, CrdS-RR468 and 
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FIG 4 Isothermal titration calorimetry analysis of HK-RR binding affinity. Thermograms and binding isotherms for the cognate CrdS-CrdA TCS pair (A), the 
noncognate CrdS-RR468 pair (B), and the noncognate HK853-CrdA pair (C) are shown. Enthalpy and Kj for binding were determined using isothermal titration 
calorimetry in a VP-ITC microcalorimeter as described in Materials and Methods. The top portion of each panel shows baseline corrected thermograms. The 
bottom portion shows the corresponding binding isotherms generated using nonlinear binding models. 
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TABLE 2 Dissociation constants determined by isothermal titration calorimetry" 

Phosphotransfer and 



Protein target 


Ligand 




phosphatase specificity 
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HK1190 


CrdA 
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HK1190 


FrgC 


>35 




HK1190 


NtrC1189 


1.23 ± 0.37 


Yes 


HK1190 


T. maritima RR468 


>55 




T. maritima HK853 


CrdA 


>75 




T. maritima HK853 


FrgC 


>55 




T. maritima HK853 


NtrC1189 


>55 




T. maritima HK853 


r. maritima RR468 


1.16 ± 0.18 


Yes 


CrdS6 


CrdA 


>75 




CrdS6 


NtrC1189 


2.09 ± 1.11 


Yes 



" Binding affinities (dissociation constants) for 4 distinct TCS HK-RR protein interactions, witli cognate and noncognate partners, were determined by ITC. Cognate TCS pairs that 
were found to interact and display phosphotransfer and phosphatase specificity are in bold. Noncognate proteins were found to have values above the upper limit detectable in 
this assay. 

^ Values are the averages of four separate measurements ( ± standard deviations}. 



HK853-CrdA, exhibited K^s above the Hmit of detection (greater 
than 75 /liM), indicating no detectable binding (Fig. 4B and C). 
Similar results were obtained for each pairwise interaction be- 
tween noncognate partners where the was above the limit of 
detection (Table 2). One caveat for this assay is that interactions 
can occur where AH is near 0 at a given temperature and would 
therefore be undetectable via ITC under the conditions tested. 
Thus, we performed additional ITC experiments at a different 
temperature (10°C). At this lower temperature, there was still no 
evidence for binding between noncognate proteins (data not 
shown). To provide further support for the above results, we mea- 
sured changes in the intrinsic fluorescence in CrdS upon binding 
to CrdA. CrdS contains a native tryptophan residue near the pre- 
dicted binding site for CrdA and thus could be utilized in fluores- 
cence binding assays. Fluorescence binding indicated that CrdS 
and CrdA interact with an apparent of 2.3 /nM, validating the 
ITC data (see Fig. S4 in the supplemental material). Together, the 
results are consistent with the hypothesis that binding affinity cor- 
responds directly to kinetic preference for both phosphotransfer 
and phosphatase function. 

Specificity residues dictate preferential binding for two- 
component systems. The results above provide strong support for 
the hypothesis that binding affinity dictates protein interactions 
for cognate signaling systems. Based on the previous work dem- 
onstrating that phosphotransfer specificity is controlled by a dis- 
crete set of residues within the HK and RR (36, 44), we wished to 
directly assess the role of those specificity residues for their effect 
on binding affinity for a subset of proteins tested above. Therefore, 
we replaced those residues within CrdS to mimic those of HKl 190 
(Fig. 5) and predicted that the chimeric protein CrdS/HK1190 
(referred to here as CrdS6) would display specificity similar to that 
of HKl 190 for its cognate partner, RR NtrC1189. We tested this 
possibility by determining the binding affinity for CrdS6 and 
NtrCl 189 via ITC and also assayed kinetic preference using phos- 
photransfer and phosphatase assays. 



To identify the correct specificity residues, we aligned the rel- 
evant sequences for CrdS and HK1190 spanning the DHp al and 
al helixes, as described previously (27). Only two of six specificity 
residues are identical between HKl 190 and CrdS (Fig. 5A). Using 
site-directed mutagenesis, we changed the remaining four speci- 
ficity residues within CrdS (resulting in the substitutions P379A, 
M382L, E385A, and T386R) to generate the CrdS6 mutant pro- 
tein. We verified that CrdS6 was active as a kinase as measured by 
the ability to autophosphorylate using ATP (see Fig. S5 in the 
supplemental material). To confirm that CrdS6 mimicked the 
phosphotransfer specificity for RR NtrC1189 rather than CrdA, 
we utilized the phosphotransfer profiling assay. CrdS6~P was in- 
cubated with equimolar amounts of each NtrC homolog as de- 
scribed above. CrdS6 was found to mimic the specificity of 
HKl 190, as was evident by a clear preference for phosphotransfer 
to NtrC 1189 without significant phosphotransfer to CrdA 
(Fig. 5B). To quantify the observed kinetic preference, we mea- 
sured the relative k^JK^ ratios as described previously (36). 
CrdS6 displays a 10^-fold increase in kinetic preference for 
NtrC1189 relative to CrdA, similar to the kinetic preference dis- 
played by wild-type CrdS for CrdA (6). In the reciprocal phospha- 
tase assay, CrdS6 was unable to dephosphorylate CrdA~P 
(Fig. 5C) yet could dephosphorylate NtrCl 189~P (Fig. 5D). Thus, 
the results clearly indicate that the chimeric CrdS6 protein 
switches its kinetic preference from CrdA to NtrCl 189. 

We then measured binding affinity using ITC to test pairwise 
interactions between CrdS6 and NtrCl 189 or CrdA (Table 2). The 
CrdS6-NtrC1189 interaction displayed a dissociation constant of 
2.1 ± 1 iM, while no measurable binding was observed between 
CrdS6 and CrdA. Moreover, the for CrdS6-NtrC1189 is close 
to the observed for each of the cognate pairs HKl 190- 
NtrC1189 and CrdS-CrdA (Table 2). Furthermore, we tested 
binding between CrdS6 and CrdA at 10°C and observed no inter- 
action (data not shown), confirming that CrdS6 lacks the capacity 
to bind CrdA. Lastly, using intrinsic tryptophan fluorescence, no 
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FIG 5 Switching specificity determinants alters binding, pliosptiotransfer, and phospliatase affinity. (A) Domain topology of^ ttie prototypical tiistidine kinase 
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binding was observed between CrdS6 and CrdA {K^ > 200 jiiM), 
further indicating these proteins do not interact (see Fig. S4 in the 
supplemental material). Together, the data lead us to conclude 
that the CrdS6 chimeric protein mimics HKl 190 for phospho- 
transfer specificity, phosphatase specificity, and binding affinity 
with RR NtrC1189. Overall, these data indicate that TCS specific- 
ity residues dictate binding affinity, which imparts kinetic prefer- 
ence and signaling fidelity. 

DISCUSSION 

Specificity residues dictate binding affinity for two-component 
signaling systems. For two-component systems, fidelity of phos- 
photransfer is predicated on specific protein-protein interactions 
between the HK and RR. Using ITC analysis, we tested four two- 
component systems and demonstrated that cognate signaling pro- 
teins had measurable binding affinities (i<C^) of approximately 
1 /i,M (Table 2) , similar to values determined previously for EnvZ- 
OmpR, CheA-CheY, CheA-CheB, and CphA-RcpA (57-59). Im- 
portantly, we observed no measurable binding for all noncognate 
protein interactions assayed by ITC. The results from the binding 
assays correlated directly with those from both phosphotransfer 
and phosphatase assays, supporting the conclusion that signaling 
fidelity depends critically on binding affinity. 



Previous work demonstrated that phosphotransfer activity is 
dependent on a set of six specificity residues (36). Based on those 
studies, we generated a chimeric protein which resulted in redi- 
rected phosphotransfer and phosphatase specificity as well as 
binding affinity between CrdS6 and NtrC1189. Furthermore, 
CrdS6 displayed neither any measurable binding affinity nor a 
kinetic preference for CrdA in the phosphotransfer and phospha- 
tase assays. Therefore, the six specificity residues described previ- 
ously are sufficient to impart both preferential binding between 
HK-RR cognate pairs and to regulate phosphorylation and de- 
phosphorylation of the correct target. In support of this view, 
recent work demonstrated that engineering noncognate HK-RR 
pairs with eukaryotic scaffolding proteins promotes phospho- 
transfer and binding both in vitro and in vivo, even in the absence 
of the known specificity residues (60). Lastly, hybrid HKs which 
contain covalently linked RFC domains typically lack specificity 
residues, presumably because those domains are spatially tethered 
(61). 

It is known, however, that cross talk can occur when either the 
cognate kinase or response regulator is absent (37, 39, 62). Thus, 
without a preferential interacting binding partner, homologous 
TCS proteins are able to engage in cross talk, as demonstrated 
during phosphotransfer profiling (53). Similarly, we observed 
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here that CrdS could phosphorylate three additional targets in 
vitro, NtrC0172, NtrC1189, and orphan NtrC3555 (see Fig. S2 in 
the supplemental material), but only at longer time points, ap- 
proaching 1 h of incubation. Importantly, CrdS and NtrC1189 
showed no detectable interaction via ITC until the specificity res- 
idues were altered to mimic those of HK1190 (Table 2). These 
results substantiate previous findings indicating that cross talk is 
enabled in vitro but is not likely relevant in vivo. Overall, we con- 
clude that binding affinity coupled with relatively fast phospho- 
transfer kinetics is the critical parameter to effectively eliminate 
unwanted cross talk within the cell and that binding affinity is 
dictated by specificity residues for cognate two-component sig- 
naling proteins. 

Signaling fidelity for TCS in M. xanthus. As a general rule, a 
given HK displays an overwhelming kinetic preference for its cog- 
nate RR. For instance, EnvZ has a 10^-fold kinetic preference for 
its cognate OmpR over the noncognate CpxR (44). The kinase 
VanS also exhibits a 10''-fold kinetic preference for VanR over the 
noncognate PhoB (63). Previous work in our laboratory demon- 
strated that CrdS has a 10^-fold kinetic preference for the cognate 
CrdA relative to the noncognate NtrC1189 (6). Thus, in vitro anal- 
ysis by phosphotransfer profiling serves as an effective way to map 
relevant signaling pathways. 

M. xanthus possesses 27 NtrC homologs, many of which have 
been demonstrated to be involved in regulation of the develop- 
mental program (Table 1). Because many mutations within this 
family of highly similar signaling proteins led to similar develop- 
mental phenotypes, we thought it possible that cross talk may 
occur or that bona fide cross-regulation may occur, as described 
for the NarXL-NarQP systems in E. coli (64). The results presented 
here, however, strongly suggest that HK-RR systems in M. xanthus 
are insulated at the molecular level, regardless of the additional 
complexity presented by their position within a hierarchy of gene 
regulation controlling development. Our results do not rule out 
cross regulation in vivo but strongly suggest that signaling fidelity 
is maintained based on the presence of a unique set of specificity 
residues for each HK associated with a given RR within the family 
of NtrC homologs. 

The system-wide (profiling) approach used in this study al- 
lowed us to assess kinetic preferences during both kinase and 
phosphatase reactions for ten distinct M. xanthus TCS, verifying 
the fidelity of CrdSA, FrgBC, HsfBA, and SasSR in vitro. The re- 
sults allow us to conclude that specificity residues dictate binding 
affinity and phosphorylation kinetics to maintain signaling fidel- 
ity in M. xanthus. It is therefore likely that complex signaling cas- 
cades regulating the developmental process in M. xanthus display 
limited cross talk and that binding affinity is the critical parameter 
for system-wide fidelity during signal transduction. 

MATERIALS AND METHODS 

Bacterial growth and DNA manipulations. For routine cloning, E. coli 
strains were grown in Luria broth with antibiotic concentrations of 
40 /J,g/ml kanamycin or 100 fig/ml ampicillin when selection was re- 
quired. All proteins were expressed using the E. coli strain BL21(DE3) 
with the isopropyl-jS-D-thiogalactopyranoside (IPTG)-inducible vector 
pET28a (Novagen) and were grown in Terrific Broth (TB). Strains are 
listed in Table SI in the supplemental material. PGR was used to amplify 
genes encoding each HK and RR by using genomic DNA isolated from 
M. xanthus strain DZ2 as a template (65). Primers for these reactions are 
listed in Table S2 in the supplemental material. Site-directed mutagenesis 



was performed using the QuikChange site-directed mutagenesis kit (Agi- 
lent), and all constructs were verified by sequencing. 

Protein purification. AH proteins used in kinase, phosphotransfer 
profiling, and phosphatase profiling assays were purified using a standard 
batch purification method, with the exception of HK0938, SasS, and 
HK3812, which were purified using a urea solubilization protocol. For 
batch purification, each strain was grown in 100 ml Terrific Broth in a 
250-ml Erlenmeyer flask at 37°C until the optical density at 600 nm 
(ODgQQ) reached 0.6 to 0.8. Protein expression was induced with 1 mM 
IPTG. After overnight growth at 20''C, cells were collected in a Beckman 
Coulter JLA- 16 250 rotor at 5,000 X g for 10 min. Cell pellets were stored 
at — 20°C until purification. For purification, cell pellets were thawed and 
suspended in 10 ml of cell lysis buffer (25 mM Tris [pH 7.6], 125 mM 
NaCl, 5 mM imidazole, 1% [vol/vol] Triton X- 100, and 0.625 g of CelLytic 
Express from Sigma-Aldrich plus Complete Mini EDTA-free protease in- 
hibitor from Roche) and incubated for 1 h before lysates were clarified by 
centrifugation at 5,000 X g. Two milliliters of Fiis-Select cobalt affinity 
resin (Sigma-Aldrich) was equilibrated in wash buffer (25 mM Tris 
[pH 7.6], 125mMNaCl, 5 mM imidazole, 1% [vol/vol] TritonX-100) and 
incubated with each lysate for 2 h at 4°C with mild agitation. His-resin was 
then pelleted by a 1 -min spin at 200 X g, and the supernatant was dis- 
carded. The resin was washed 3 times with 1 0 ml of wash buffer and eluted 
with 5 ml of elution buffer (25 mM Tris [pH 7.6] , 125 mM NaCl, 500 mM 
imidazole). Samples were dialyzed overnight at 4''C against dialysis buffer 
(25 mM Tris [pH 7.6], 125 mM NaCl, 1 mM dithiothreitol [DTT], 1% 
[vol/vol] Triton X-100, 50% [vol/vol] glycerol, 0.5 mM EDTA). Purified 
protein was assayed by standard denaturing gel electrophoresis, and the 
concentration was determined using the Bradford reagent. AU purified 
proteins were stored in dialysis buffer at — 20°C. 

Proteins used for ITC analysis were prepared from cells grown to an 
ODjoo of 0.8 in 2 liters of Terrific Broth, induced with 0.5 mM IPTG and 
shaken overnight at 20°C. Cells were pelleted by centrifugation and stored 
at — 20°G until purification. For purification, cell pellets were thawed and 
suspended in 25 ml of ITC ceUlysis buffer (25 mM Tris [pH 7.6], 125 mM 
NaCl, 0.1% J vol/vol] Triton X-100, 1 mg/ml lysozyme. Complete Mini 
EDTA-free protease inhibitor JRoche]) and lysed by passage through an 
LV02 Microfluidizer (Microfluidics). Sample lysates were then clarified by 
centrifugation at 40,000 X g and passage through a 0.45-/j,m filter disk. 
Lysateswere loaded onto a 5-ml Hi-Trap HP (GE) column using an AKTA 
fast protein liquid chromatography (FPLC) system (GE). The column was 
washed with 7.5 column volumes (CV) of buffer A (25 mM Tris JpH 7.6] , 
125 mM NaCl, 10 mM imidazole, 0.1% [vol/vol] Triton X-100) and 
eluted with a 15-CV Unear gradient to 100% buffer B (25 mM Tris 
[pH 7.6], 125 mM NaCl, 500 mM imidazole, 0.1% Jvol/vol] Triton 
X-100). Fractions containing protein of interest were combined and dia- 
lyzed overnight against ITC buffer (25 mM Tris JpH 7.6], 150 mM NaCl, 
10% glycerol, 0.1% Jvol/vol] Triton X-100, 0.5 mM EDTA). Proteins were 
then assayed for purity using SDS-PAGE. If there were any apparent con- 
taminating proteins, samples were further purified using a Superdex S75 
(GE) gel filtration column equilibrated in ITC buffer. 

Kinase assays. Kinase assays were performed as described previously 
(16). Briefly, 5 p,l of 50 /nM HK stock was added to 5 /aI 1 0 X kinase buffer 
(250 mM Tris JpH 7.6], 500 mM KCl, 10 mM MgClj, 10 mM MnClj, 
10 mM CaClj and 10 mM dithiothreitol) and 35 p,l distilled H2O. Reac- 
tion mixtures were incubated at room temperature and started by the 
addition of 5 jml of an ATP mix (250 ijlM ATP, 3 jjlM Jy-^^PJATP). Ali- 
quots were removed at various time points and stopped by the addition of 
an equal volume of 2X SDS loading buffer. Samples were resolved by 
electrophoresis on 12% SDS-polyacrylamide gels. The dye front, contain- 
ing unincorporated ATP, was removed, and then gels were exposed for 2 
to 4 h on a phosphor screen and visualized using a Typhoon Imager (GE). 
ImageQuant v 5.1 was used to determine integrated pixel density and 
subtract background. 

Phosphotransfer profiling. Histidine kinases were allowed to auto- 
phosphorylate, as detailed above, until maximal phosphorylation levels 
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were reached (10 min to 2 h, depending on the kinase). A 5-fj.\ aliquot was 
mixed with 5 of purified RR (10 fxM stock) in kinase buffer, and reac- 
tions were stopped at various times by addition of 2 X SDS loading buffer. 
Incorporation of labeled phosphoryl groups was analyzed as detailed 
above. Each phosphotransfer profiling experiment was run on three sep- 
arate gels. Gel images were handled identically and assembled into figures 
using Microsoft PowerPoint. 

Phosphatase profiUng. For phosphatase profiling experiments, each 
RR was radiolabeled with acetyl phosphate ([^^P]AcP) as described pre- 
viously (6, 16). Each RR~P was then incubated in kinase buffer for 5 min 
with equimolar amounts of each HK (final concentration of 5 fxM for each 
HK and RR). RR~P incubated with buffer alone was set to 100%. Samples 
were separated by SDS-PAGE and visualized as detailed above. Bar graphs 
show averages and standard deviations from three experiments. 

Binding measurements using isothermal titration calorimetry. 
Measurements were performed using a VP-ITC titration calorimeter 
from Microcal, Inc. Each histidine kinase and response regulator was 
dialyzed overnight against the same ITC buffer. The VP-ITC cell was filled 
with the histidine kinase protein at concentrations of~10tol5 fiM in ITC 
buffer. The syringe was filled with -150 to 200 /aM response regulator in 
ITC buffer. The exact concentration of each protein was measured using a 
bicinchoninic acid protein assay kit (Thermo Scientific). The ITC cell and 
syringe solutions were extensively degassed. The chamber was kept under 
constant stirring at 350 rpm, and all experiments were performed at 25°C 
except where noted. Control injections were also performed to determine 
heats of dilution. The control injection data were subtracted from the raw 
data using the Origin ITC analysis package provided with the microcalo- 
rimeter. The heat of dilution not eliminated by controls was determined 
by the average of the last 3 to 5 injections and subsequently subtracted. 
The data were analyzed using the single-site binding model provided with 
the ITC analysis package. Results in Table 2 are the averages from at least 
four separate runs from two separate protein preparations. 

Intrinsic tryptophan fluorescence measurements of binding affini- 
ties. Intrinsic tryptophan fluorescence was used to monitor the binding of 
cognate and noncognate response regulators with histidine kinases. For 
each pair, the response regulator (-100 fxM) was titrated into 1.3 ml of 
5 /J,M histidine kinase contained in a stirred 1-cm quartz cuvette until 
saturation. Fluorescence intensity measurements were recorded at 25°C 
on a Fluorolog 3 (Jobin Yvon, Fioriba) spectrofluorimeter. For each mea- 
surement, excitation was at 294 nm (5-nm slit width). The fluorescence 
emission spectrum was collected in the range 300 to 400 nm (5-nm slit 
width). Individual measurements were integrated for 30 s. An equal vol- 
ume of buffer was titrated into the histidine kinase, and intensity mea- 
surements were collected under the conditions described above. These 
measurements were then subtracted from the corresponding response 
regulator titration to correct for the change in intensity due to dilution. 
The for each titration was determined by fitting buffer-corrected data 
using nonlinear regression (Sigma Plot; SPPS, Inc.). The resulting disso- 
ciation constant was the average of three independent experiments. 
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